Calculation of the P-T phase diagram and tendency toward decomposition in 

equiatomic TiZr alloy 
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Electronic, structural and thermodynamic properties of the equiatomic alloy TiZr are calculated 
within the electron density functional theory and the Debye-Griineisen model. The calculated 
values of the lattice parameters a and c/a agree well with the experimental data for the a, to and ft 
phases. The uj phase is shown to be stable at atmospheric pressure and low temperatures; it remains 
energetically preferable up to T = 600 K. The a phase of the TiZr alloy becomes stable in the range 
600 K< T < 900 K, and the ft phase at temperatures above 900 K. The constructed phase diagram 
qualitatively agrees with the experimental data available. The tendency toward decomposition in 
the equiatomic alloy ui— TiZr is studied. It is shown that in the ground state the ui phase of the 
ordered equiatomic alloy TiZr has a tendency toward ordering, rather than decomposition. 

PACS numbers: 63.20.Ry, 05.10.Gg, 63.20.Kr, 71.15.Nc 
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Introduction 

It has been experimentally found that the TiZr sys- 
tem is characterized by full solubility of its constituents. 
As in pure titanium and zirconium, three phases (a, (3 
and uj) are observed in the TiZr alloji~— . The structural 
a — > (3 transformations of the equiatomic alloy were ex- 
tensively studied in Ref. 1 by the differential thermal anal- 
ysis (DTA) at temperatures up to 1023 K, and pressures 
up to 7 GPa. It was found that the f3 — > a transition tem- 
perature, being equal to 852 K at atmospheric pressure, 
decreases with pressure down to the triple equilibrium 
point of the a, (3 and w phases (P tr = 4.9 ± 0.3 GPa, 
Ttr — 733 ±30 K). At pressures above the triple point 
the f3 phase transforms immediately to the id phase with 
a light positive slope of the equilibrium line. If a sample 
is cooled to room temperature at a pressure of 6 GPa, 
and then unloaded, one can obtain at atmospheric pres- 
sure a metastable id phase which on heating transforms 
into an a phase in the temperature interval from 698 K 
to 743 K. Cooling of the (3 phase in the pressure range 2.8 
- 4.8 GPa results in the formation of a two-phase mixture 
of a stable a and a metastable id phase. The structural 
a — > lo transformations in the TiZr alloy were studied in 
detail in RefA Investigating TiZr samples under shear- 
strain conditions at pressures up to 9 GPa at tempera- 
tures 300 and 77 K, the authors arrived at the conclusion 
that in equiatomic TiZr the equilibrium a — > id boundary 
is situated on the P-T diagram at 6.6 GPa. In the same 
paper the phase diagram of TiZr was constructed in the 
regular-solution approximation, and the triple point pa- 
rameters were calculated (P = 8.5 GPa, T = 693 K). As 
may be seen, these values differ substantially from those 
of Refi 

Detailed studies performed in Ref. 4 have shown that 
in the region of high pressures and temperatures there 
exist two id phases (u) and u>\) that differ in atomic vol- 
ume by about 14%. The authors suggested the existence 
of a isostructural phase transition co — u\ connected with 



changes in the electron structure of the alloy. They sup- 
posed that the large difference in the atomic volume be- 
tween the two phases points to the existence of an s-d 
electronic transition in to— TiZr. Later, phase separation 
of a hexagonal TiZr w phase was experimentally detected 
in RefA The id — > id% + cl>2 decomposition was revealed 
after a prolonged heat treatment at P — 5.5±0.6 GPa 
and T = 440 ± 30° C. It was supposed that in a wide 
concentration range at pressures above the triple equilib- 
rium point, the id phase may exist in the Ti x Zr\- x alloy 
only as a metastable one that persists due to low diffu- 
sive mobility of its constituents. The decomposition of 
the id— TiZr solid solution into two id phases of different 
structure was used as an alternative explanation for the 
experimental results obtained in Ref. 4 . 

Up to now the electron structure and structural trans- 
formations of the equiatomic alloy TiZr have never been 
calculated. Below we present the results of our theoret- 
ical calculations of electronic, structural and thermody- 
namic properties of the equiatomic alloy TiZr performed 
within the framework of the electron density functional 
theory and the Debye-Griineisen model. The tendency 
of the ordered equiatomic alloy u— TiZr to decompose is 
also investigated. 



I. CALCULATION TECHNIQUE 

The electron structure and total energy were calcu- 
lated by the scalar relativistic full-potential linearized 
augmented-plane-wave (FPLAPW) method, using the 
WIEN2K package 7 . To ensure the desired accuracy of 
the total energy calculation, the number of plane waves 
was defined by the condition RK max — 7, the total num- 
ber of /c-points in the Brillouin zone was equal to 3000, 
3000, 600 for the /?, a and lu phase, respectively. The 
total and partial densities of states were obtained by a 
modified tetrahedron method^. The atomic radii were 
the same for all phases and pressures: 2.42 a.u. for Zr, 
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and 2.26 a.u. for Ti. 

In FigQ]are shown the crystal structures of the j3, a, 
and oj phases of the equiatomic alloy TiZr used in the 
calculation. It is seen that the /3 phase was represented 
by a structure of the CsCl type with Ti atoms at the cube 
sites and a Zr atom at the center. The a phase had a 
hexagonal close-packed lattice in which one atom was Zr, 
the other Ti. Finally, to describe the oj phase we used a 
hexagonal lattice with 6 atoms per cell (an oj structure 
doubled along the z axis). The atomic arrangement and 
species in this case were chosen as follows: (0,0,0) - Zr, 
(3' §> 4) " Ti, (3, 31 4) - Zr, (0, 0, 5) - Ti, (3, |, |) - Zr, 
(|, i, |) - Ti. Thus, the j3 and a phases were represented 
by layers of Ti and Zr alternating along the z axis, and 
in the oj phase layers of Ti and Zr were separated by a 
mixed Ti-Zr layer. 
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FIG. 2: Volume dependence of the total energy of equiatomic 
TiZr alloy for the f3, a and oj phases 
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FIG. 1: Crystal structure of the f3, a and oj phases of 
equiatomicTiZr. The Zr atoms are grey, the Ti atoms are 
dark. 



In the hexagonal structures the ratio c/a was optimized 
for the experimental volume values. In the following, 
when calculating the volume dependence of the total en- 
ergy, the c/a ratio was considered to be constant. In 
Table [TJ are listed the calculated and experimental equi- 
librium values of the TiZr lattice parameters. 



TABLE I: Equilibrium values of the TiZr lattice parameters 
in atomic units 
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1.583 
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OJ 


9.122 


0.617 


9.152 


0.617 


9.131 


0.616 



It is seen from the table that the lattice parameters a 
and c/a obtained in our calculation agree well with the 
experimental data. The greatest discrepancy is observed 
for the lattice constant of the oj phase. It should be noted 
that in Ref— the lattice parameters were calculated for a 
metastable oj structure at atmospheric pressure, and in 
Ref4 for pressure-strained samples. 

The total energy of each structure was calculated for 7 
values of the cell volume V. The data obtained were then 
interpolated using the technique proposed by Moruzzi 8 . 
Such an interpolation scheme, together with the Debye- 
Griineisen model, makes it possible to include implicitly 



anharmonic effects. The curves obtained for the volume 
dependence of the electron subsystem total energy are 
shown in Fig. [21 The energy zero in figure corresponds 
to -8906.0 Ry. 

As seen from the figure, the energy minimum in the 
ground state falls on the oj phase. And only at the 
relative volume change V/Vq — 0.75 the j3 phase be- 
comes energetically preferable. A similar situation was 
observed in pure Ti and Zr as well. However the dif- 
ference in energy between the a and oj phases in TiZr 
amounts to A£" Q _ W = 6 mRy, while in pure Ti and Zr it 
is 0.8 mRy and 1 mRy, respectively. Hence it follows that 
in the equiatomic alloy TiZr the stability region of the oj 
phase should be much larger in temperature than in pure 
metals Ti,Zr. The equilibrium values of the volume are 
V eq = 20.47 A 3 , 20.52 A 3 , and 19.97 A 3 for the a, oj and 
P phases, respectively. 

Figure [3] presents the volume dependence of the free 
energy at different temperatures. The free energy of the 
TiZr alloy was calculated in the Debye-Griineisen model 
with allowance made for the contributions from the elec- 
tron entropy. The technique of calculating this latter has 
been described in detail in Ref<2. As seen from the figure, 
the relationship between the energies of different struc- 
tures changes with temperature. So, at 300 K and zero 
pressure the energy minimum falls, as in the ground state, 
on the oj phase which remains energetically preferable up 
to T = 600 K. In the temperature range 600 K<T<900 K 
it is the a phase of TiZr alloy which becomes stable, while 
above 900 K it is the f3 phase. 



II. PHASE DIAGRAM 

The phase diagram of TiZr based on the analysis of 
Gibbs potentials for different structures is presented in 
FigHl The results of calculation are shown by the solid 
line. The dotted line denotes the experimental equilib- 
rium boundaries for the a, (3 and oj phases of TiZr ob- 
tained in Ref. 1 . The experimental values of the a — oj 
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FIG. 3: Free energy of the a, lo and /3 phases of TiZr at 
different temperatures 



transition at room temperature are taken from Refs^i 3 -. 

On the whole, a good agreement of the calculated triple 
point {Ptheor — 4.2 GPa, T t heor = 720 K) with the exper- 
imental values P exp = 4.9 ± 0.3 GPa, T exp = 733 ± 30 Ki 
is observed. At zero pressure the calculated temperature 

943 K. This value is 



of the P - a transition is Tf*° r 
higher than the experimental one, Tt x ^ a — 852 K, defined 
in Refi as the average of the temperatures of the tran- 
sition onset on heating and cooling. It should be noted 
that a large hysteresis is observed upon the a — f3 trans- 
formation in TiZr. At atmospheric pressure the maxima 
of thermal peaks in the DTA curves fall on T ~ 912 K 
on heating and T ~ 810 K on cooling, the typical peak 
width being AT ~ 40 K. With this in mind, one can 
consider the results of calculation of the a — (3 equilib- 
rium boundary in the Debye-Griineisen model as quite 
satisfactory. 

The greatest discrepancy between the theoretical cal- 
culation and the experimental evidence available is ob- 
served for the a — uj transition. In Ref£ the pressure 
at which this transition occurs at room temperature was 
estimated to be P^J p bJ = 6.6 GPa. Note thet equilib- 



rium point of the a and uj phases was determined under 
shear-strain conditions at pressures up to 9 GPa. The 
shearing strain is known to lower the pressure at which 
the phase transition begins. Presumably for this reason 
the authors of Ref. 3 failed to obtain the a — oj transition 
at room temperature under quasi-hydrostatic conditions. 
In Refj^ it was shown by X-ray diffraction method that 
the a phase of TiZr remains the sole stable phase under 
quasi-hydrostatic pressure up to 12.2 GPa. Only from 5.5 
GPa on, becomes dominating the uj phase which remains 
stable up to 56.9 GPa. At pressures above 56.9 GPa there 
forms a high-pressure phase with a bcc lattice. 
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FIG. 4: The P-T phase diagram of TiZr. The solid line shows 
the calculation results. The dotted lines are constructed from 
the experimental data [l|. The experimental values for the 
a — uj transition at room temperature are taken from Refs.: 

• a. 

As seen from FigUl in our calculation at atmospheric 
pressure and low temperatures the uj phase is stable, 
there occurs noa-w transition at room temperature. 
Note that in our calculations of pure Tii£ and Zr— , in 
complete agreement with the experimental data, the a 
phase is stable at atmospheric pressure and room tem- 
perature, and the uj phase is stable only under pressure. 
That the uj phase in TiZr at normal conditions is energet- 
ically preferable, immediately follows from a comparison 
of the calculated free energies (see FigJ3]). Recall that 
the difference in energy between the a and uj structures 
in the equiatomic TiZr alloy is almost five times greater 
than in pure titanium and zirconium. 

The discrepancy between experiment and theory may 
be due to the fact that the calculation was performed for 
ideal crystalline structures (see Fig[T]) , whereas the exper- 
imental samples were imperfect crystals with lattice de- 
fects. In particular, it was shown 3 - that in a TiZr alloy 
shear-strained under pressure the uj phase is represented 
by aggregations of oblong particles with characteristic 
size of 3—5 nm, and 15—30 nm long. If uj— phase particles 
are situated in a coarse grain of a phase, they are mainly 
located at its boundaries. It was also noted^ that various 
imperfect structures in samples pre-treated in different 
ways have a noticeable effect on the course of structural 
transformations. Evidently, we could not model a real 
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structure in first-principles calculations. 

The correctness of our results may be supported by the 
following experimental evidence^: firstly, the metastable 
lu phase was obtained at atmospheric pressure as a result 
of cooling of the (3 phase under a pressure of 6 GPa with 
subsequent unloading at room temperature. Secondly, 
between 2.2 and 4.8 GPa on cooling of the f3 phase there 
forms a two-phase mixture of a stable a and a metastable 
lu phase. And lastly, it was found that at atmospheric 
pressure the u phase in the TiZr alloy, when heated above 
698 K, transforms into an a phased. In Ref.- the tem- 
perature of this transformation was defined as T = 623 
K for P = 0.0001 GPa. This value differs by only 13 K 
from the temperature Tp=Q = 610 K we have calculated 
for the uj — > a transition. 




FIG. 5: The structure types used in modeling the decompo- 
sition in uj — TiZr. The Ti atoms are dark, the Zr atoms light. 



II two layers of pure titanium alternated along the z axis 
with two layers of pure zirconium with no intermediate 
layer. For structure III were chosen six Ti layers alternat- 
ing with six Zr layers without intermediate layer. And 
lastly, in variant IV five Ti layers were separated by an 
intermediate layer from five Zr layers. 

The free energy was calculated by the scalar rel- 
ativistic full-potential linearized augmcnted-plane-wave 
(FPLAPW) method, using the WIEN2K package!. In 
the first two variants the number of atoms per unit cell 
was six (3 Ti atoms and 3 Zr atoms). In variants III 
and IV the number of atoms in the cell amounted to 18 
(9 atoms of each species). For variants I and II struc- 
tural optimization of the c/a ratio was performed, and 
the equilibrium atom positions were defined with the pro- 
cedure of minimizing the forces acting on atoms. 
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The above evidence suggests, in our opinion, that in 
the TiZr alloy with ideal crystal lattice the phase dia- 
gram should look as it is depicted in FigJU Of course, 
it must be taken into account that the Debye-Griineisen 
model, used for calculating the thermodynamic poten- 
tials, is a rather rough approximation and, obviously, 
cannot ensure good accuracy, especially at high temper- 
atures when the anharmonic effects become of consider- 
able importance. It should be noted that at room tem- 
perature the pressure calculated for the uj — > (3 transition 
is nearly half as large as the experimental value. Because 
of low temperatures, we do not believe this discrepancy 
to be connected with the choice of the Debye model for 
describing the thermodynamic properties. It is rather 
due to the deviation of real alloys from the ideal peri- 
odic structures used in calculating the total energy in 
the ground state. 

III. CALCULATION OF THE TENDENCY 
TOWARD DECOMPOSITION 

To estimate the tendency toward decomposition in 
uj— TiZr in the ground state, the total energy was cal- 
culated for four structures (see Fig [5]). Structure I was 
represented by layers of pure titanium and pure zirco- 
nium alternating along the z axis and separated by in- 
termediate mixed Ti-Zr layers. In the two-layer structure 



FIG. 6: Energy change upon lattice relaxation with respect 
to the equilibrium atom positions in ideal uj phase for variants 
I and II 



In Figj6] the volume dependence of the lattice relax- 
ation energy change AE — E re i ax — Eq is displayed for 
calculation variants I and II. Here Eq is the total energy 
of the system with atomic arrangement corresponding to 
the ideal uj lattice; E re i ax is the system energy after min- 
imization of the forces acting on atoms for a given vol- 
ume. As seen from the figure, with decreasing volume the 
atoms become displaced from the positions correspond- 
ing to the ideal uj lattice, the displacement magnitude 
depending on the volume and the structure type. For 
the structure of type I, corresponding to the most uni- 
form distribution of Ti and Zr atoms at V < 15 A 3 , there 
occurs a sharp decrease of AE due to significant atomic 
rearrangement. For the two-layer system (II) such a re- 
arrangement is not observed in the considered interval of 
volume change. In FigLZl f° r V — 13 A 3 are depicted 
the (110) planes, the arrows indicating the direction of 
atomic displacements on relaxation for the lattices of type 
I (a) and II (b) (the z axis is pointing upwards) . It is seen 
that in both cases the atomic displacements are directed 
only along the z axis. In the two-layer system (II) the 
atoms of titanium and zirconium are displaced in op- 
posite directions, whereas in system I the atomic chain 
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displacement occurs without strain. The displacements 
shown in Fig(7]^a) correspond to those characteristic of 
the 10 f3 transition. The volume value V « 15 A 3 at 
which begins a sharp decrease in AE, agrees well with 
the results of the total energy calculation for the (O and 
/3 phases of TiZr plotted in FigEJ Based on the data 
presented, we can draw an important conclusion that 
the pressure value at which occurs the u> — > (3 transi- 
tion depends substantially on the ordering type in the 
equiatomic TiZr alloy. This also indirectly confirms our 
statement that the disagreement with the experiment 
concerning the to — /3 equilibrium boundary position on 
the phase diagram calculated in the Debye model (Figj4]) 
is caused by the presence of inhomogeneities in actual 
TiZr alloys used in experiments in Refsji~— . 



(a) (b) 




FIG. 7: The direction of atomic displacements upon lattice 
relaxation for calculation variants I (a) and II (b) 



The volume dependence of the total energy for the re- 
laxed structures of type I and II is plotted in FigEJa). 
As may be seen, in the whole range of volume change 
the energy is minimum for the structure of type I corre- 
sponding to the most uniform distribution of Ti and Zr 
atoms. Recall that in variant I the lattice is represented 
by pure monolayers of titanium and zirconium separated 
by a mixed Ti-Zr layer, and in variant II by a system 
of two Ti layers alternating with two Zr layers (with no 
intermediate layer). Thus, as decomposition grows, the 
system energy increases. This tendency persists on fur- 
ther decomposition, which may be seen in FigJSJb), where 
the total energy is plotted versus volume for the struc- 
tures of type II, III and IV. In the structure of type IV 
five Ti layers are separated from five Zr layers by a mixed 
intermediate layer, while in structure III there are six lay- 
ers of each metal without intermediate layer. It is seen 
from the figure that as the thickness of pure metal layers 
increases (from one to five layers) , the system energy sig- 
nificantly rises. As could be expected, the presence of an 
intermediate layer reduces the total energy of the system. 
This may be seen from a comparison of the energy values 
at V ~ 18.5 A 3 for the calculation variants III and IV 
(structure III is presented in the figure by a single point). 
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FIG. 8: Total energy of ui— TiZr for different types of decom- 
position ( I - IV) designated in accordance with Fig[5] 



To summarize, the calculations performed show that in 
the ground state the uj phase of TiZr exhibits a ten- 
dency toward ordering and not toward decomposition, as 
was suggested in Ref.—. The analysis of the total energy 
curves shows that the allowance for temperature effects 
in the Debye-Griineisen model will not change the en- 
ergy relation between different calculation variants, and 
cannot explain the experimentally observed formation of 
two uj structures. Besides, we have not found any pe- 
culiarities connected with the s-d electron transition in 
the total-energy curves. Thus, also the suggestion ad- 
vanced in Ref. 4 that there exists an isostructural transi- 
tion lo — > lo\ due to the pressure-induced changes in the 
electron structure, is not confirmed by the calculation. In 
our opinion, the high-temperature decomposition in the 
uj phase of equiatomic TiZr alloy is connected not with 
the change in electron structure under pressure, but with 
peculiarities of the lattice dynamics, in particular, with 
the presence of strongly anharmonic vibrational modes 
which are of crucial importance in stabilization of the uj 
lattice of pure titanium and zirconium^. 
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